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ABSTRACT 


The time dependent mixed-layer depth and temperature re- 
sponse of a two-layer hydrothermodynamic ocean model to a 
moving hurricane model was investigated. The hurricane model 
was that of Elsberry, Pearson, and Corgnati (1974). The 
‘maximum eye radius and hurricane intensity was bounded to 
isolate ocean reaction from hurricane variance. The ocean 
model was based on the upwelling model of O'Brien and Hurl- 
burt (1972) in two dimensions and was located along the path 
of the storm. The momentum and depth equations in this 
model are treated semi-implicitly. Mechanical and convective 
mixing was included and was found to be the dominate cause 
of deepening of the mixed layer and cooling of the ocean sur- 
face. Alternating regions of upwelling and downwelling were 
produced in the wake of the storm as a result of inertial 
waves caused by the passage of the storm. The cases studied 
indicated a linear relation between the wavelength and the 


storm translation speed. 
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I. INTRODUCTION 


It is generally accepted that sea-surface temperatures 
may be strongly affected by tropical storm activity, and 
vice versa. Sea-surface temperatures in excess of 26°C are 
apparently required to produce hurricane-force winds. It 
has also been observed that cooling of as much as 6°C may be 
present in the wake of a hurricane. 

One of the first systematic efforts o investigate the 
cooling effect was done after the passage of hurricane Hilda, 
1964. From this data, Leipper (1967) found clear evidence 
supporting previous scattered observations which indicated 
lowering of sea-surface temperatures in the wake of a hurri- 
cane. The sea-surface temperature was observed to drop 5 to 
6°C from the temperature present before the storm passage. 
The observations indicated upwelling from about 60 m in depth. 
Associated with this was an outward transport of water which 
occurred in the warm surface layers with cooling and mixing 
of these layers as they moved outward. Finally this outward 
moving water experienced convergence and downwelling around 
the hurricane area. The observations were taken along paths 
perpendicular to the path followed by the storm. This would 
not allow certain features along the path of the storm to be 
observed, for example, the internal waves predicted by 


Geisler (1970). 
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In conjunction with project Stormfury, Black and Mallin- 
ger (1972) reported on data taken before, during, and after 
the passage of hurricane Ginger, 1971. These oceanic mea- 
surements were taken from a Navy WP-3 aircraft using air- 
borne expendable bathythermographs, (AXBTs). The primary 
information gained, of interest in this thesis, was that as 
the storm moved faster less surface cooling occurred. In 
fact, there appeared to be a cut-off speed related to the 
ratio of the baroclinic wave speed and the speed of propaga- 
tion of the storm. Geisler (1970) predicted this ratio to 
be Y2 whereas the data obtained suggested a ratio of 3. At 
the slower speeds strong upwelling near the center should 
occur. While, at the faster speeds the upwelling intensity 
decreases, and therefore does not contribute greatly to the 
surface cooling. 

Sheets (1974) presented evidence indicating internal 
gravity waves at the thermocline in the wake of a rapidly 
moving hurricane, Ellen, 1973. A total of 18 AXBTs were 
dropped at intervals of five to fifteen nautical miles. In 
addition to the internal wave indication, a surface cooling 
of .5°C to 1.5°C was observed. The amount of surface cool- 
ing is in agreement with the concept that a rapidly moving 
storm does not cool the ocean surface to the same extent as 
slower moving storms. In the same paper, however, data from 
hurricane Irah, a slow moving storm (7-10 kt), showed nomi- 
nal cooling, less than .5°C, instead of the greater antici- 


pated cooling. These data illustrate that many factors play 
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A IA ne response Ol the ocean surface temperature to 
the passage of a hurricane. 

O'Brien and Reid (1967) developed a model of the time- 
dependent response of a two-layer (one active layer) ocean 
model under a Stationary, axially-symmetric atmospheric 
vortex. The dynamic response of the ocean was also assumed 
to be axially-symmetric, and therefore solutions in two in- 
dependent variables, radial distance and time, could be ob- 
tained. Gradient winds were obtained from a pressure 
distribution function resembling actual data. The surface 
stress was derived using these winds. The results of this 
model indicated strong upwelling near the center due to ve- 
locity divergence. This upwelling was strongest near the 
region of maximum surface stress. After 48 hours the depth 
of the upper layer had been reduced to about 10 meters in 
the region of maximum upwelling. 

Geisler (1970) developed a model of the linear response 
of a two-layer ocean to a stress pattern representing a 
hurricane moving at a constant speed. When the speed of 
movement was greater than the baroclinic wave speed, this 
model produced internal waves at the interface of the mixed 
layer and the thermocline in the wake of the storm. The 
wavelength of these disturbances increased as the speed of 
movement of the storm increased. The disturbance propagated 
in the same direction and at the same speed as the storm. 
Another feature of this model was that a baroclinic ridge 


was left in the ocean in the wake of the storm. Due to 
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upwelling near the center and downwelling in the surrounding 
area, a wake of reduced mixed-layer depth near the center of 
the path was produced. The amplitude of this ridge was 
greater with decreasing speed of movement of the hurricane. 

The models discussed thus far treated only the response 
of the ocean to a given wind stress. In an effort to devel- 
op a model to describe the interaction of the ocean and the 
atmosphere during a hurricane, Elsberry, Pearson and Corgnati 
(1974) developed an axisymmetric hurricane model. As the 
atmospheric portion of this model is used in this thesis it 
will be described later. The ocean used in this case was at 
a fixed temperature. 

Fraim (1973) introduced a mixed-layer model capable of 
producing time-dependent values of sea-surface temperature 
and mixed-layer depth in response to atmospheric forcing. 

It was based on the mixed-layer models of Kraus and Turner 
(1967), and Denman (1973), which predicted the mixed-layer 
temperature and depth in response to wind mixing and convec- 
tion. The radial current was defined by the Ekman current 
response to the stress at the surface. In this model the 


stress input to wind mixing, T was a function of the total 


0 > 
stress, Ta? the roughness length, Z,, the mixed-layer depth, 
ho, and an adjustable constant C. 


ho 
CZ.) 





Toug Ta EXP (- 


The excess stress, T. - T was used to produce currents from 


a 0? 


Ekman balance considerations. This yielded a depth-averaged 
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lilii current, and due to symmetry, this is the only cur- 
rent which could produce upwelling. Application of the 
continuity equation to these values of current gave the up- 
welling at the bottom of the mixed layer. This model was 
used to study the effects of the three cooling processes in- 
volved in lowering the sea surface temperatures. These are, 
entrainment mixing across a stable layer, convective mixing, 
and vertical and horizontal advection. The results of these 
experiments showed that both advection and mixing have sig- 
nificant effects on the sea surface temperature and mixed 
layer depths, and as expected these effects are most pro- 
nounced in the region of maximum surface stress. 

Trapnell (1974) introduced more refinements into the 
Fraim model, the most notable of which was the addition of a 
return flow in the lower layer and the ability to let the 
Surface stress move across the ocean. The flow in the lower 
layer was assumed to have the same mass transport as the 
upper layer. This added a new method by which the below 
layer gradient could be changed. In addition to upwelling 
and downwelling, there now existed horizontal advection. 
Movement of the stress and cooling functions was simulated 
by translating the center of the atmosphere grid relative to 
the ocean grid. This adjustment was made at certain time 
increments to simulate different speeds. The results from 
Trapnell's experiments suggest that entrainment mixing is 
the primary mechanisms for cooling the oceans, except within 


the region of maximum winds where intense upwelling becomes 





important. Also it is seen that these effects are strongly 
associated with slower moving storms. Due to the lack of 
dynamics this model did not exhibit the internal waves pre- 
dicted by Geisler. 

One of the primary constraints in both Fraim's and 
Trapnell's models was that the currents used were derived 
from Ekman balance to a given surface stress. As a further 
refinement, this thesis describes a two-layer hydrothermo- 
dynamic model and its response to surface stress and ae 


representing a moving hurricane. 
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JI. MODEL 


A. OCEAN MODEL 

The dynamical equations in the ocean model and the meth- 
od of solution described in this thesis follow the coastal 
upwelling model of O'Brien and Hurlburt (1972). However the 
present model contains a term to include thermal and mechan- 
ical mixing due to upward heat flux and mechanical energy 


input at the surface. 
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Equations (1) and (2) are the radial momentum equations 
in the upper and lower layers respectively. These include 
the advection, pressure gradient, coriolis, stress, and 


lateral diffusion terms. 

















S I 
OV} OV} = T Yr y das, 
St + Du fu, + oh, + ATZ2 (3) 
Ly By 2 
IV> O ns: T -T + IV) A 
at i ar ph» ne > 


17 





These are the tangential momentum equations, comprised of the 
same terms as the radial momentum equations, excepting the 
pressure gradient forces, which due to symmetry are zero. 


The stress in the momentum equations is defined by 
E = -k Ta sin a (5a) 


TY = k T, COS Q (5b) 


where k is the partitioning constant which determines how 
much of the surface stress will be used for current genera- 
tion, and a is the atmospheric inflow angle. Following 
O'Brien and Hurlburt (1972) the stress components of the 


interface are defined as 


I 
1% =p cq(ui-u2) (6a) - 
I e 
1% =p cq(vi-v2) (6b) 
where q = (qitq2)/2 
neon 
and q; = (u;+v;) 


The constant c is the drag coefficient at the interface, and 


is chosen to be 3x10 * also following O'Brien and Hurlburt. 
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Equations (7) and (8), which describe the temperature and 
depth of the mixed layer, are the same as used by Fraim 
(1973) and Trapnell (1974). The term M is the contribution 
to depth change due to entrainment mixing. The increase in 
potential energy of the column through entrainment at the 
interface is equal to the excess of mechanical generation 
(G) of kinetic energy over dissipation (D), plus a contribu- 
tion due to convective overturning 1f the column is losing 
heat (ptg > 0). A is called the entrainment parameter 
which is 1 when M is positive and zero when M is negative. 
Thus mixing across the base of the mixed layer occurs only 
when M is positive, and therefore acting to deepen and cool 
the mixed layer. If M were permitted to be negative it would 
Simulate heat extraction from the lower layer, to warm the 
already warmer upper layer. As this model depicts a wind 


dominated region, M is always i. 











dh. , Shou, _ _ 

ot + ðX M > 
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Equations (9) and (10) are the equivalent equations to 
(3) and (4) but in the lower layer. Mixing does not contrib- 
ute directly to the change in temperature at that level, 
therefore only the advection and vertical motion terms are 
present. It can be shown that with the assumed linear dis- 


tribution of currents and temperatures in the lower layer, 
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that the temperature advected by the average current compo- 
nent iS that from a level 2/3 of the distance from the bot- 
tom of the lower layer. 


Ju» 
ðX 





W; = —ho (11) 


T. = 2T, - T 


h (12) 


b 


Equations (11) and (12) are diagnostic equations used to 
close the equation set. Equation (11) was derived by inte- 
grating the continuity equation in the lower layer and using 
the boundary conditions of no mass flux through the bottom 
of that layer. Equation (12) utilizes the fixed bottom tem- 
perature, the linear lapse rate assumption, and the newly 
computed values of the depth and average temperature of the 
lower layer to compute a new value of temperature at the top 
of the thermocline, Tp: 

The boundary conditions allowed no heat or momentum flux 
through the ends or bottom of the model. At the surface a 
momentum flux from the atmosphere to the ocean was applied. 
A heat flux from the ocean to the atmosphere also occurred 
at the surface. The model was considered to be a free sur- 
face model because the height of the ocean surface varied 
with respect to the fixed bottom depth. 

As in O'Brien and Hurlburt the solutions to the equa- 
tions contain external and internal gravity waves. There- 
fore the pressure gradient terms in the radial momentum 


equations were treated implicitly to allow a longer time 
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step to be taken and still retain these waves. Treatment of 
these equations is shown in Appendix A. The remainder of 
the equations were evaluated using centered differencing in 
time and space. Stability problems will occur due to using 
centered differences in both time and space in the diffusion 
term, therefore time lag was used in this term. 

The atmospheric parameters were computed at 100 grid 
points with a three kilometer spacing. This atmospheric 
forcing was moved at a constant speed across this grid in 
order to simulate the motion of a storm. The “em phe 
forcing was applied to the ocean model as shown in figure 1. 
Between and including C and D the atmospheric forcing param- 
eters were applied directly, from D to E the values at D 
were applied. The values from C to D in the atmosphere were 
applied to the points from C to B in the ocean, although the 
sign was changed on the stress terms for symmetry. From 
point B to point A the values at point B were applied. 
Points B, C, and D move from right to left to simulate the 
motion of the storm. 

As in Trapnell's model the surface stress and heating 
functions were not initially applied at full value. To per- 
mit the initially quiet ocean to adjust, the forcing was 
increased exponentially, requiring about nine hours for ap- 
plication of 95% of the full force of the storm. It was 
found that both the atmospheric and oceanic models remained 


stable for time steps of 900 seconds. 
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B. ATMOSPHERIC MODEL 

The atmospheric hurricane model used to drive the ocean 
in these experiments is the two-layer, time-dependent model 
of Elsberry, Pearson, and Corgnati (1974) which was also 
used in the moving model by Trapnell (1974). The hurricane 
was considered to be mature so that an axially-symmetric 
storm could be assumed. The maximum winds were permitted to 
intensify to 50 m sec”! and then held fixed. This allowed 
the changes in the ocean structure to be isolated from 
changes in the atmospheric model. 

The atmospheric model was integrated simultaneously dur- 
ing these experiments using the new sea-surface temperatures 
computed by the ocean model. New surface stresses and sur- 
face heating values were also applied to the ocean. This 
constituted the interaction between the ocean and the atmos- 
phere. 

This model consisted of a surface inflow layer and an 
upper outflow layer. The tangential wind flow at the surface 


outside the radius of maximum winds was given by 


ble 


V r< = conse 


and inside the radius of maximum winds by 


Var’ 2 const. 


These equations were empirically derived by Riehl (1963) and 
supported by observations by Gray and Shea (1973). The 


radius of maximum wind was assumed to coincide with the 


23 





eyewall radius. Normally in this model the eyewall radius 
is defined to be that radius at which the average of five 
adjacent values of equivalent potential temperature was a 
maximum. However, to isolate the changes in the ocean, the 
eyewall radius was fixed at 30 km after a short period of 
adjustment. The boundary layer in the atmospheric model was 
a two-layer, baroclinic model developed by Cardone (1969). 
The heat and momentum fluxes at the interface were computed 
in this model. These were a function of the stability of 
the boundary layer and the wind at the top of this layer. 
To obtain more information regarding the atmospheric model 
the paper by Elsberry, Pearson, and Corgnati (1974) should 


be consulted. 
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III. RESULTS 


A. EFFECTS OF TRANSLATION SPEED 


The simulated storm was moved at various speeds, and 


compared with the stationary case to observe the changes in 


the behavior of the ocean due to storm translation Speed. 


Many of the model results are consistent with previous 


theories and observations, as discussed earlier. The be- 


havior of selected ocean parameters is displayed in figures 


3 through 6 and 10 through 13. The parameters displayed in 


these figures are: 


a) 
b) 


c) 


d) 


€) 


f) 


the temperature at the top of the thermocline; 

the temperature of the mixed layer, and therefore 
the ocean surface: 

the depth of the mixed layer; 

the vertical velocity as derived from the continuity 
equation in the lower layer; 

the upper-level current component directed along the 
path of the storm (positive values indicate a cur- 
rent toward the right; also note that the baseline 
IS 1 m sec); 

the tangential wind velocity at the top of the atmos- 
pheric boundary layer (the stress outside the storm 
is held at the value present at the outer edge of 
the storm, rather than being zero as implied by 


this graph). 
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In the standard model half of the surface stress was 
used to produce currents, and half was used for energy input 
to entrainment mixing. The total surface stress had a maxi- 
mum value of 2.6 nt m”? (26 dynes cm”?). The mixed layer 
had an initial depth of 30 m and temperature of 30°C, while 
the lower layer, the thermocline, extended 100 m below the 
mixed layer and terminated with a fixed bottom temperature 
of 20°C. This initial state is depicted in figure 2. As 
can be seen from figure 2, a small ledge in the temperature 
appears at the bottom of the mixed layer. From Eq. (7) it 


is seen that if To is equal to T, the term M becomes infi- 


h 
nite. To maintain a slightly stable layer at the interface, 
this difference is initialized at .1 degree and never al- 
lowed to become smaller in magnitude. However, as this tem- 
perature difference immediately increases in response to the 
wind mixing, the condition was not forced following initial- 
ization. 
1. Stationary Storm Model 

The results after integration for 36 hours with the 
translation speed equal to zero are displayed in figure 3. 
The wind profile is shown in figure 3f. Note that this does 
not depict the wind outside the storm area. The vertical 
motion at the bottom of the mixed layer is displayed in fig- 
ure 3d. Combining this with the radial currents in the 
upper layer, as shown in figure 3e, we see there are two 
circulation cells on each side of the storm center. The 


first cell has its maximum upwelling at the inside edge of 
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the radius of maximum winds, RMW, and downwelling at the 
outer edge of the RMW. A second cell of about half the 
magnitude is located at about 90 km from the center. The 
maximum depth of the mixed layer is observed at the outer 
edge of the RMW, rather than coincident with it as might be 
expected. At the RMW the deepening due to mixing is offset 
by upwelling. 

Referring to figures 3b and 3d the minimum tempera- 
ture of the mixed layer does not correspond to the maximum 
upwelling region, or from figure 3c, to the maximum depth of 
the mixed layer. Instead the lowest temperatures are di- 
rectly under the RMW, indicating that both upwelling and 
mixing contribute to cooling of the upper layer. The tem- 
perature at the top of the thermocline is shown in figure 
3a. Deepening of the mixed layer, and the resultant reduc- 
tion in thickness of the lower layer, would be the greatest 
factor in the lower temperature at the top of the lower iay- 
er. The reduced thickness would lower the temperature di- 
rectly due to the lower position on the thermal lapse rate. 
The water being advected into this region would also be 
colder, thus contributing to lower temperatures (see the 
discussion of Tj, in the model equation (10)). 

Table I compares Trapnell's Ekman current model with 
this dynamic model. The results are for a stationary storm 
case after being integrated 18 hours. The major difference 
is in the minimum mixed layer depth. This could be due to 


the slower reaction of the dynamic model, as this value was 
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about 24 m after 27 hours. Another explanation could be the 
different method of partitioning the surface stress. From 
this comparison the results from a model using currents in 
Ekman balance appear to compare favorably with the dynamic 
model. 
2. Moving Storm Model 

The results of moving the storm at 3.2 kt for 72 
hours are shown in figure 4. Comparing figure 4b with fig- 
ure 3b, it is seen that a minimum mixed-layer temperature of 
28.3°C was observed as compared to 27.3°C associated with 
the stationary storm. Evidence of waves on the interface at 
the bottom of the mixed layer is shown in figure 4c. The 
wavelength is about 180 km. These waves are more noticeable 
in the P oal and upper layer radial velocities in figures 
4d and 4e. 

The comparison of Trapnell's Ekman model moving at 
3.39 kt and the dynamic model moving at 3.2 kt is shown in 
Table II. The main difference is the momentum terms, how- 
ever overall these values of maximums and minimums do not 
differ greatly. The problem with the Ekman model is that it 
does not contain the dynamics necessary to produce the hori- 
zontal distribution of variables that are present in figure 
4. It appears that for sea-surface temperature prediction 
the Ekman balance model would produce comparable results to 
the more complicated dynamic model. However for a predic- 
tion of the internal structure of the ocean, (e.g. mixed- 
layer depth, slope of the thermocline, etc.) the dynamic 
model is necessary. 
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The translation speed of the storm was increased 
from 3.2 kt to 6.5 kt. In figure 5 the results of inte- 
erating this model 36 hours is presented. This can be com- 
pared with figure 4 as the same distance has been traveled 
by the storm but in one-half the time. Note the longer 
wavelength in the faster storm and also the lag of the up- 
welling region behind the storm center is much longer in 
this case. Notice also, comparing 4b and 5b, that less 
cooling has occurred due to the passage of the faster storm 
with a minimum mixed-layer temperature of 28.9°C compared 
with 28.3°C in the slower storm. At the top of the thermo- 
cline the same results occur with 26.5°C for the faster 
storm and 25.0°C for the slower storm. This case was also 
integrated to 72 hours and these results are presented in 
figure 6. This will be considered the standard case for the 
following discussion. Notice that there is little differ- 
ence compared with the 36 hour case in figure 5, except for 
the longer wake behind the storm. Temperatures were only 
slightly lower, by .3°C for the top of the thermocline and 
-1°C for the mixed layer. The wavelength and magnitude of 
the variations in mixed-layer depth, vertical velocity, 
radial velocity and temperature at the top of the thermo- 
cline remain essentially the same. This wavelength is ap- 
proximately 360 km as compared to 180 km in the slower storm. 

These results may be compared with the data from 
Sheets (1974). In the observations the temperature at the 


interface is higher for greater mixed-layer depths, MLD, and 


36 





ERMOCL INE 


ee 


| 
| 


T 


O 
VS 


MIXED LAYER 


Ra 
AR 


O 
no 


Man Cayce 


a 


(€) 
2 


O) 
(J 


(0) 


e 


(SD 
Gi 


(M) 


ar 


Figure 5. 


SI ads 


3UU 600 + 900 1200 L200 
Per a) 


36 Hour Predicted Values With the Storm 
Moving at 6.35 kt. 





a 


OCEAN VERTICAL 


aa 
D 


OCERN HADIAL 
VEL. 
= 


=, 
Nu? 


Een n 


[r (o) 


Chie 
dez 


ur) 
Oy 


een 


UY 


í 


s. NY usa 


A ES 


500 SOU 


DISTANCE UH] 


Figure 5. 36 Hour Predicted Values With the Storm 
(Cont.) Moving at 6.5 kt. 


38 





THERAMOCL INE 


IA 


o 
Nx 


Mao A 
mee 


O 
< 


A DAA 


=a 


CD 
OS 


aca 


Lid 


Si 


ee) 
2 


Oy 


ME 
A 


J 


Figure 6. 


STARMI CENTER 
301 600 900 1200 1500 


IS) 


72 Hour Predicted Values With the Storm 
Moving at 6.5 kt. 


39 





US 


` = 


po 


SL co 
(SPBA-N3 “130 
D aaa ND 


3 


m~ 


QO 


ae 


00 


| 
Es 
= 
= 
(3 
pe 
3 _ 
5 
D 
Oo 
Ss 
cr e 
x”. 
|- 
= 
Ly 
A A A, Á. 
=> 
y en C3 
(Gis cr) 
E 
— 
un 
A Oy oD 
A~ Dee Sé oO 


MESA 2 
Widbs hiii 


(J4AS7W5 GNI™M 
oH JBIINJSNOL 


Bone a. 


72 Hour Predicted Values With the Storm 


Moving at 6.5 kt- 


Figure 6. 


(Cont. ) 


40 





lower for shallow MLD's, which implies that the waves are 
due to upwelling and downwelling rather than mixing. In 
figure 7, derived from the data of the standard case we see 
that this same relationship occurs. However, in Sheets! 
data these temperature fluctuations also appear within the 
mixed layer, whereas these fluctuations are not present in 
figure 6b. Also in the observations the wavelength is on 
the order of 35 to 40 km which does not agree with the model 
wavelength of about 360 km. Geisler's theory predicts the 
longer wavelengths, about 185 km for a 5.2 knot storm. More 
data with different storm speeds are necessary to compare 
with the theoretical and model predictions of wavelengths. 
The vector contributions to acceleration of currents 
in both the upper and lower layers are presented in figure 
8. These values are computed at the center of the region. 
The sum of the vectors is proportional to the instantaneous 
acceleration at that point. The first set of vectors, at 27 
hours, represents the ocean condition prior to passage of 
the storm. At this time the forces are very nearly in bal- 
ance, except for the lower-level current component normal to 
the storm track (v2). The storm passes this point after 37 
hours. However, at 36 hours the stresses have been reduced 
to a very small value compared with the pressure gradient 
and coriolis terms. It is seen that this sudden reduction, 
and subsequent reversal of the stress pattern, produces 
strong inertial waves in the currents. It should also be 


noted that the phase of the pressure gradient term also 
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contributed to the imbalance of accelerations. At 45 hours 
this term and the coriolis term are in the same direction, 
thus contributing to a large imbalance. This condition 
would not occur without the lag of the pressure gradient 
term behind the center of the storm. After 45 hours it is 
seen that the coriolis acceleration dominates the solutions, 
exchanging kinetic energy between the current components 
directed parallel (u) and normal (v) to the storm track. To 
better illustrate this, a hodograph at point 250 of the cur- 
rent direction and magnitude as a function of time is shown 
in figure 9a. Before the storm passes, the winds are from 
the top of the figure producing a stress which creates cur- 
rents from that direction, and due to coriolis there is a 
current away from the storm center. This is the same behav- 
ior noted in the stationary case. In the lower layer the 
pressure gradient term produces currents toward the center 
with the coriolis acceleration producing a current component 
perpendicular to the path of the storm. The storm passes 
from the right at about 37 hours. The inertial effects are 
then seen to begin, with the current direction reversed be- 
tween 36 and 45 hours. Note that the inertial period for 
the position of the storm, 25°N, is 28.34 hours, and observe 
that from t=45 hours to t=72 hours, almost one rotation of 
the current direction has occurred in 27 hours. Therefore 
we see that inertial effects are very strong, and indeed 
continue to dominate after the passage of the storm center. 


weto this relation to the inertial period, and because the 
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waves move with the same speed and direction as the storm, 
it becomes apparent that the wavelength should be greater 
for faster moving storms. This relationship is linear in 
the cases studied here, that is when the storm movement 
speed is doubled the wavelength is doubled. 

A hodograph for the lower layer at this same point 
and for the same storm speed is presented in figure 9b. The 
currents in the lower layer have the same features as the 
upper layer except the magnitudes are smaller and are 180° 
out of phase with the upper layer. Remembering that the 
vertical motion is a function of the first derivative of the 
currents in the lower layer, it is seen that vertical cur- 
rents are 90° out of phase with the fluctuations of lower 
layer current parallel to the storm track. This upwelling 
and downwelling produces the waves at the bottom of the 
mixed layer, therefore their wavelength is also the same as 


the fluctuations in currents. 


B. VARIATION OF OCEAN PARAMETERS 

To check the sensitivity of the model, the standard run 
was repeated while varying individual parameters. In this 
section the effects of increased stability in the thermo- 
cline layer, increased initial mixed-layer depth, increased 
stress for currents, and the elimination of the mixing term 
will be discussed. 

Changing the temperature of the undisturbed layer, and 


hence the temperature at the bottom of the thermocline, 
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changes the stability of that layer. In this experiment, 
shown in figure 10, the temperature was lowered to 15°C, 
thus increasing the stability. The major differences were a 
slightly cooler mixed layer, about .1°C cooler than the 
standard model, and a significantly cooler temperature at 
the top of the thermocline, 1.8°C cooler than the standard 
case. The reason for the relatively small change in the 
mixed-layer temperature was that it required more energy to 
lift the cooler subsurface water into the mixed layer through 
mixing. The significant reduction in the temperature at the 
top of the thermocline is due to the cooler water being ad- 
vected into the region due to the steeper lapse rate (refer 
to equation (10)). The currents and depths remained essen- 
tially the same. 

The next variation from the standard case was to ini- 
tialize the mixed-layer depth at 50 m. These results are 
Shown in figure 11. Due to the larger amount of mass pres- 
ent within the layer due to the increased mixed-layer depth, 
the velocities in the upper layer were smaller, with a maxi- 


! compared to .48 m sec-! for the standard 


mum Of .36 m sec” 
case. Also, the mixed layer was .5°C warmer in this case. 
Again the wave structure remained unchanged although the 
overall depth of the mixed layer was increased. 

The results of changing the partitioning of surface 
stress to 70% for current production and 30% for mixing are 


Shown in figure 12. AS expected, the currents were stronger 


Petits case than in the standard case of 50% partitioning 
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shown in figure 6. Comparing figure 6b with 12b it is seen 
that the surface temperature is lower when the mixing is 
stronger. Figures 6a and 12a show the same to be true for 
the temperature at the top of the thermocline. This result 
and the deeper mixed layer in 6c compared to 12c, show that 
mixing is the most important contributor to deepening of the 
mixed layer. A close comparison of the waves at the bottom 
of the mixed layer, however, shows that the amplitude is 
increased by about 6 m. It therefore appears that the am- 
plitude of these waves is a function of the magnitude of the 
currents. 

A run was also made with the mixing term set equal to 
zero. Comparing figure 13 with those of the standard case, 
figure 6, it is seen that in the case without mixing the 
surface and thermocline temperatures have only been slightly 
affected. The mixed layer has not deepened significantly 
although the waves were present with about the same ampli- 
tude. The currents were also stronger even though the 
stress component partitioned to current production is the 
same. This effect is due to the mixed layer remaining 
Shallow, and therefore the magnitude of the depth-averaged 


currents are larger. 
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IV. CONCLUSIONS 


The results shown in the previous section agree quite 
well with previous theories and observations regarding the 
behavior of the ocean beneath a moving hurricane. The cool- 
ing of the surface temperature, as observed by Leipper 
(1967), Black and Mallinger (1972), and Sheets (1974), was 
produced by this model. Entrainment mixing was the dominate 
factor in cooling of the surface temperature and deepening 
of the mixed layer. Upwelling maxima and minima as pre- 
dicted by Geisler (1970) were also produced in the wake of 
the moving storm. Increased translation speed produced less 
cooling behind the storm; this agreed with the data of Black 
and Mallinger. In comparison with Trapnell (1974) this 
model appears to be quite consistent thermodynamically. 
However, as stated previously the earlier model did not con- 
tain the dynamics to produce the horizontal structure pro- 
duced by the present model. 

The primary discrepancy was the wavelengths of the 
waves at the bottom of the mixed layer. Sheets shows values 
on the order of 50 km for a fast moving storm whereas this 
model shows wavelengths on the order of 360 km for storms 
moving at 6.5 kt. This difference has not been resolved. 

In this model these waves appear as a result of inertial 
waves produced by the sudden change in stress as the storm 
passes and due to the lag of the pressure gradient force 


behind the storm center. 
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In addition to resolving the difference in wavelength 
mentioned above, other work is necessary on this model. The 
method of partitioning stress for mixing and current produc- 
tion needs to be studied. Turbulent momentum transfer at 
the interface should be included. Sensitivity tests should 
be run on the drag coefficient at the interface and the dif- 
fusion constant. 

The most valuable rerinenent would be to extend this 
model to three dimensions. This would eliminate the most 
restrictive assumption in this model, axial symmetry. To 
make this modification, and maintain the present resolution 
in the most active region, near the hurricane center, would 
require some type of nested grid arrangement. The nested 
erid approach would be necessary due to the present limita- 


tion of computers in speed and size. 





APPENDIX A 
SEMI-IMPLICIT SCHEME FOR 


MOMENTUM AND DEPTH EQUATIONS 


The equations in finite difference form which were used 
in the computer model are developed here. Using the equa- 
tions presented in the discussion of the model we proceed to 
met solutions, for u, u,, h, and h,. Following the 


procedure in O'Brien and Hurlburt (1972) 
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The tendencies are now evaluated for a 2At time step and the 
divergence terms are averaged between time levels (n-1)At 
and (n+1)A4t. The subscript j is the space index while the 


Superscript n is the time index. 
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n+1 ð n+1 
uj E 





= “E (A-5) 
a atm, er =D Th 2 ZN er (A-6) 
“a ee) IX j 
upy +lAtggy(bitha)-Atg' ELIT = zatu, +u," 
-[Ates c(h tha )-atg SUIT” = La OD 


petty an 


== = 
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n n-1 
= ar -[AtH25 L; (A-8) 


With homogeneous boundary conditions on u, and u2, h, and h; 
can be eliminated. Substituting (A-6) and (A-8) into (A-5) 
and (A-7) respectively, and dropping the ntl superscript 


the following is obtained. 





d gu; 9U2 
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‚+{AtgS-[L2-(AtH, Sl 


T. )+Ly-(AtH2 S44) ) 








-atg Zu [12-(Attızal)]}, = (A-10) 


Expanding (A-9) gives: 





late AH] = Li-Atggzlhotlu) (A-11) 
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Expanding (A-10) gives: 


dU} = 








AGA MI (A) GA 





ð t 
= La-Atg[2(Lz+L,)]+Atg 52 (A-12) 





Subtracting (A-11) from (A-12) gives: 
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MU te (UNE Hs ) = L3-L,+Atg (A-13) 
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Substituting uij from (A-13) into (A-11) gives: 
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For the remainder of this derivation only, the following 


definitions are made. 


a = gHı (At)? 
b= gH: (At)? 
© = g Hi(At)? 
CAS)" 


From (A-13) we get 


3 u 





u oe = e (A-15) 
From (A-14) we get 
92 
era asu, = _ 
U; (a DE a d (A 16) 
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eee rl 
Me °. 4441 =i 
Using ae = EE Eq. (A-15) becomes 


e m = 7 
u; j y E e 24, 40, 3-1? e (A-17) 


E. Al " a eee A-18 
Uy 54,7 Ct SU; ü,j-1 == ( ) 


This is the form to apply the method presented in Richtmyer 


(1967). Using the symbol conventions in that text, 


A = 1 

E rl: 
b = (2 + cf? 
C= 1 
D 


uns Te de 


As c and f are both positive the criterion, explained in 
Richtmyer, of B>A+C is satisfied. 
Having used a guess for Uy; above, a new a 3 using the 


above Ur; can now be computed from Eq. (A-10). 


U, 11621 rr = = i ee u) (A-193) 
As used for u, 

A = 1 

B=2+ 

C = 1 

D = + A, 41720, 40,31) 
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Solve Eq. (A-18) and (A-19) iteratively to obtain values for 
u, and uz. 


The scheme in Richtmyer is 


u. = E.u. ,+F. A-20 
j I ae 
where 
A. 1 
B.-C.E. B.-E 
i J J l! J j= 
and 
D F. 1 
tae 7 (A-22) 
J B.-E. 
eee 
The boundary conditions, ea eee Oommen, are Zero. 


Values for Ong are obtained by rearranging Eq. (A-6) 


n+1 = 3 dU, oa a 
na Lo [AtHi 5 E (A-23) 





Likewise Eq. (A-8) provides a relationship for Dart 


nri du». ntl 
2 j Ly- [Athos 1; (A-24) 





h 


This completes the semi-implicit scheme used in computing 
new values for u, u>, h, and h>. 
The remainder of the equations are handled explicitly, 


using centered differences in time and space. 
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